ID to 2D Na Ion Diffusion Inherently Linked to Structural Transitions in Nao.7Co02 
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We report the observation of a stepwise "melting" of the low-temperature Na-vacancy order in 
the layered transition metal oxide Nao.7Co02. High-resolution neutron powder diffraction indicates 
the existence of two first-order structural transitions, one at Ti ^ 290 K, followed by a second 
at T2 ~ 400 K. Detailed analysis reveals that both transitions are linked to changes in the Na 
mobility. Our data are consistent with a two-step disappearance of Na-vacancy order through the 
successive opening of first quasi- ID (Ti > T > T2) and then 2D (T > T2) Na diffusion paths. These 
results shed new light on previous, seemingly incompatible, experimental interpretations regarding 
the relationship between Na-vacancy order and Na dynamics in this material. They also represent 
an important step towards the tuning of physical properties and the design of tailored functional 
materials through an improved control and understanding of ionic diffusion. 

PACS numbers: 71.27.+a, 71.30.+h, 71.45.Lr, 61.05.fm 



Amongst correlated electron systems, layered transi- 
tion metal oxides (TMO) have attracted a particular in- 
terest in condensed matter physics for their intriguing 
magnetic [1] and electronic [2| properties. The key to 
their understanding is believed to lie in the fundamen- 
tal physics of the TMO planes, with the structural ion- 
units sandwiched between them acting merely as passive 
charge reservoirs. In parallel, many of these materials 
are in the center of attention for applied sciences, in par- 
ticular in the field of rechargeable batteries [3|. Here 
the focus is instead on the intermediary ionic layers and 
in particular their dynamic properties, while the role of 
TMO planes is considered as secondary. Very recently 
these two fields have been unified under the framework of 
the layered Naa^Co02 compounds 0]. While its magnetic 
[HI and electronic [6] properties change dramatically with 
the number of holes/Co (i.e. x) in the Co-0 layers 0, 
it has lately become evident that also thepotential land- 
scape created by ordered Na vacancies [1, Q as well as the 
Na dynamics [lO, UjJ play equally important roles. This 
has hinted at the possibility of tuning physical properties 
through the control of the self-assembling ion-texture 
(i.e. periodic potential landscape) on a sub-nanoscale. 
The key-issue is here to understand the solid state dy- 
namics at the atomic level. 

Below we present a high-resolution neutron powder 
diffraction (NPD) study which clearly reveals how Na 
dynamics and changes in the crystal structure are inher- 
ently linked in the Naa;Co02 system. For this purpose 
we have chosen Nao.7Co02, whose Na-vacancy ordered 
array is amongst the most stable of the NaajCo02 family 
[m . The symmetry of this pattern has been investigated 
in the past but the actual periodicity as well as 

the size and arrangement of Na-vacancy clusters are still 



subject of debate. The temperature dependence of this 
array has been much less investigated, but a few groups 
have predicted [15] as well as reported 
existence of a first order transition at Ti 



ill, HI 



UJ, 116[ the 
290 K. The 

nature of the high temperature phase (T > Ti) is so far 
contradictory. Single-crystal neutron diffraction studies 
suggest the stabilization of a static stripe- like ordered Na 
vacancy array [l3l|. In contrast, NMR JjJ and /i+SR in- 
vestigations [l0| indicate an increasing Na diffusion for 
T > 200 K and a "melting' of the low-temperature Na- 
vacancy array above Ti. 

With the aim of clarifying these two apparently incom- 
patible results we have re-investigated the crystal struc- 
ture of Nao.7Co02 as a function of temperature. Approx- 
imately 10 grams of powder sample were synthesized as 
described in the Supplemental Material [SM] . A series of 
NPD patterns (T = 10 - 475 K) were recorded at the 
high resolution diffractometer HRPT [13] (monochroma- 
tor: Ge (533), A = 1.494 A) of the Swiss Neutron Spal- 
lation Source SINQ, Paul Scherrer Institut, PSI Villigen, 
Switzerland. Scattering from the sample environment 
(cryo-furnace) was suppressed using an oscillating radial 
collimator. All data were analyzed using the Rietveld 
refinement package of the Full Prof Suite [18]. 

The NaajCo02 family is characterized by a layered 
structure with alternating C0O2 slabs and interstitial, 
partially occupied, Na^^ layers [SM]. Two crystallo- 
graphic, non-equivalent Na sites are located at the nodes 
of two interpenetrating triangular lattices: Nal, which 
lies above and below the Co sites in the neighboring C0O2 
layers, and Na2, located at the center of a Co trigonal 
prism. The average crystal structure of Nao.7Co02 has 
been described using a primitive unit cell of P6s/mmc 
symmetry (a/, ^ 2.83 A, Ch ^ 10.88 A [SM]). Ex- 
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tra superstructure reflections have been observed above 
and below Ti and attributed to different Na-vacancy 
ordered patterns. The corresponding super cehs have 
been reported to be commensurate with the underlying 
P63 / mmc lattice which, according to [sl, ITU, , remains 
metrically hexagonal for T = 1.5 — 350 K. 

The present data reveal new details concerning the 
metrics of the average unit cell. The temperature de- 
pendence of the high-angle region in the neutron diffrac- 
tion patterns [Fig.IlJa)] clearly shows the presence of two 
flrst-order structural phase transitions involving changes 
in the PGs/mmc fundamental reflections. The flrst of 
them occurs at Ti ^ 290 K, coinciding with the anomaly 
described in previous studies [III, [HI, [l6| . A second tran- 
sition, not previously reported, is found at higher tem- 
peratures (T2 ^ 400 K). The neutron diffraction pat- 
terns above T2 can indeed be reflned using the P63 / mmc 
space group [SM] . This is however not possible for the in- 
termediate region between Ti and T2, where a splitting 
of the hexagonal reflections with nonzero Miller indices h 
and/or k is clearly visible [see Fig. [TJa)]. This indicates 
that the a.h and hexagonal axis are inequivalent within 
this temperature range. The diffraction patterns can be 
indexed on an orthorhombic cell 3.0 = s.^, ho = hh - 
Co = c/^ (a/i and hh denote here pseudohexagonal axes) 
and reflned using the Cmcm space group, which is one 
of the maximal non-isomorphic subgroups of PGa/mmc. 
The relationship between the hexagonal and orthorhom- 
bic cells is schematically shown in the inset of Fig. [IJb). 
Tables with the atomic positions are given in [SM]. 

The splitting of the a.^ and hh orthorhombic axes at 
T2 upon cooling gives rise to an expansion of the ab- 
plane, which coincides with a much weaker increase of c 
[SM]. As shown in Fig.[T]^b) the splitting decreases again 
below Ti, although it remains nonzero at the lowest in- 
vestigated temperature. The Rietveld flts using either 
the P6s/mmc or Cmcm space groups are signiflcantly 
worse in this temperature region. In particular, the extra 
splitting/broadening of some of the high- angle reflections 
cannot be reproduced without decreasing the symmetry 
from orthorhombic to monoclinic [SM]. The fits using 
the P2i/m space group and a unit cell very similar to 
that of Cmcm with 7 slightly larger than 90° result in a 
clear improvement. However, the agreement between the 
observed and calculated intensities is worse than in the 
intermediate and high temperature regions (see differ- 
ence patterns and values in [SM]). This indicates the 
existence of subtle structural features which can not be 
reproduced by the average P2i/m monoclinic cell. Such 
behavior is consistent with the existence of complex Na- 
vacancy ordering below Ti. 

We shall show now that the occurrence of this "quasi- 
reentrant" structural rearrangement lies in the connec- 
tion between structural changes and Na-ion dynamics in 
Nao.7Co02. The evolution of the in-plane distances be- 
tween neighboring Na ions is shown in Fig.[2fa). The low- 




Temperature [K] 

FIG. 1: (Color online) (a) 2D contour plot showing the tem- 
perature dependence of the neutron powder diffraction pat- 
terns for Nao.7Co02. (b) Relationship between the in-plane 
hexagonal (h), orthorhombic (o) and monoclinic (m) lattice 
parameters together with its thermal evolution. The vertical 
lines indicate the two transition temperatures Ti ^ 290 K 
and T2 ~ 400 K. N.B. The error-bars are smaller than the 
marker-size. 

temperature monoclinic structure used to refine our data 
allows three distinct, although nearly indistinguishable 
Nal-Na2 distances I (long), m (medium) and s (short), 
see the insets of Fig.[2](a) and [SM]. At Ti two of them be- 
come suddenly identical (5) and shorter (-2.5%) whereas 
the third one {I) undergoes a dramatic increase (+5.2%). 
The difference between I and s changes very little be- 
tween Ti and T2. Then, it drops to zero at the transi- 
tion to the high temperature hexagonal phase (T > T2), 
where the three Nal-Na2 distances become identical (m). 

The changes of the Nal-Na2 distances at Ti and T2 are 
about one order of magnitude larger than those under- 
gone by any other interatomic distances in the structure. 
This can be appreciated in Fig. EJa), where the tem- 
perature dependence of the three Co-0 distances is also 
displayed. Particularly noticeable is the absence of signif- 
icant discontinuities in the c lattice parameter since, ac- 
cording to other studies, displacements of the Co and/or 
O atoms along this crystallographic axis would be ex- 
pected as response to the changes in the Nal/Na2 occu- 
pations. This lack of anomalies is not due to compensa- 
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FIG. 2: (Color online) Temperature dependence of (a) the Nal-Na2 and Co-0 distances, (b) the thickness of the 0-Na-O layer 
and (c) that of the C0O2 slabs, (d-f) Nal and Na2 site occupations. Inset in (a) is a schematic representation of the Na-triangle 
distortions and Nal-Na2 lengths (l, m, s) for the different space groups. Solid lines and shadings are guides to the eye. 



tion effects, as demonstrated by the smooth temperature 
dependence of the thickness of the C0O2 slabs and the 
0-Na-O layers [Figs. [2jb-c)]. We thus conclude that the 
main coherent atomic displacements occurring at Ti and 
T2 are confined to the Na-ion planes. 

The existence of a correlation between structural 
changes and Na dynamics is first indicated by the ther- 
mal evolution of the Nal/Na2 fractional occupations, ni 
and n2 which, together with the total Na content ni + 
712 is shown in Fig. [2fd-e). At T = 10 K, the ener- 
getically more favorable Na2 position displays an occu- 
pancy approximately twice that of Nal, in agreement 
with previous reports [l9|. The refined value of ni + 
712 0.7) coincides well with the x value derived from 
the lattice parameter c [lH and it does not change be- 
low Ti within experimental error. This is not the case 
for Til and n2, whose temperature dependence indicate a 
weak but progressive Na transfer from the Na2 into the 
Nal sites above ~T=200K. At Ti a clear anomaly is ob- 
served that coincides with a substantial increase of the 
Na2 Nal transfer. Such behavior is consistent with 
the growing diffusion above T = 200 K and the "melt- 
ing" of the Na.y layers at Ti suggested by NMR [ll| and 



/i+SR studies [lO[. An interesting detail is that the sum 
of the Na content at the Nal and Na2 sites (ni + TI2) 
remains nearly constant below Ti . However, it under- 
goes a small but significant decrease (~ —0.02) at Ti 
[Fig. El^e)]. This means that, for T > Ti and within the 
interaction time of thermal neutrons, approximately 3% 
of all the Na ions are located away from the Nal and/or 
Na2 positions [Fig. Eff)]. 



Further experimental evidence suggesting an active 
role of Na dynamics is provided by the temperature de- 
pendence of the Debye- Waller (DW) factors. These quan- 
tities reflect the reduction of Bragg integrated intensities 
due to the displacements of atoms from their equilibrium 
positions. For a harmonic crystal, they can be inter- 
preted in terms of time-averaged mean-square displace- 
ments <u^> resulting from all normal modes of vibra- 
tion. However, the <u^> values obtained from structural 
refinements often contain information about other types 
of displacements, either static or dynamic, which are not 
included in the structural model. 

The evolution of the mean-square displacements un 
and U22 at the Nal site along the a and b axes, respec- 
tively, (monoclinic/orthorhombic notation) is displayed 
in Fig. [3](a). At T = 10 K U22 is nearly zero, consistent 
with the reduced thermal motion expected at this tem- 
perature. Around T = 200 K, however, a rapid increase 
suddenly sets on. In contrast, un is rather large already 
at T = 10 K and displays only a weak T— dependence, 
signaling a static contribution superimposed onto a very 
weak thermal motion. This is consistent with the exis- 
tence of Na-vacancy order, proposed in previous studies 
and suggests a distribution of Nal positions along the 
direction. For the Na2 site, large values of un and U22 
below Ti along with a weak T— dependence indicate once 
again the presence of static displacements but, contrary 
to Nal, the associated distribution of Na2 positions is 
nearly isotropic (i.e. un ^ 1^22)- 

At Ti clear anomalies are observed in the Na mean- 
square displacements un and 1^22- As shown in Fig. [3l 
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FIG. 3: (Color online) Temperature dependence of the in- 
plane mean square displacements un and U22 for the Nal and 
Na2 sites, as well as the out-of-plane U33 for the O site. Top 
insets display the shape of the projection of the associated 
Nal/Na2 thermal ellipsoids in the ab-plane 



the in-plane projection of the Nal and Na2 thermal el- 
lipsoids drastically change their shape. This may arise 
either from a new distribution of static displacements, 
e.g.^ from a new Na- vacancy arrangement [l3| or. as 
suss; ested by NMR [uj and /i+SR measurements |1Q| . 
from the onset of Na diffusion. Based on the DW 
factors alone it is not possible to distinguish between 
these two possibilities. However, the changes in the 
Nal /Na2 fractional occupations mentioned above clearly 
favor the second one. A further arguments is that the 
largest projection of the Nal/Na2 thermal ellipsoids is 
now observed along the shortest Nal-Na2 distances {s) 
whereas the contribution along / is significantly smaller 
[Fig. IIJd]. The mobility of diffusing ions in ionic con- 
ductors is known to largely depend on the connectivity 
between the available hopping sites. Hence, this obser- 
vation suggests the activation of quasi-one-dimensional 
(QID) ...Nal^^^Na2^^^Nal^^^Na2... zigzag dif- 
fusion pathways at Ti coinciding with the shortening of 
2/3 of the Nal-Na2 distances [Fig.[2K and Fig.©)]. Above 
T2 a similar correlation between Na+ motion and struc- 
tural changes can be established: the thermal ellipsoids 
of all Na sites become perfectly isotropic coinciding with 
the equalization of Nal-Na2 distances (m) . This sug- 
gests a crossover from QID to 2D Na+ diffusion at T2. 

Further experimental evidence supporting this scenario 
is given by the Fourier difference maps [SM] which pro- 
vide information about the residual scattering density 
contributing to the Bragg reflections not reproduced by 
neither the structural model. Three representative exam- 
ples of the results obtained for the Na-ion layer {z = 0.25) 
at T = 50, 350 and 450 K are shown in Fig. 11 with the 
Nal-Na2 framework superimposed (see also [SM]). At low 
temperatures [Fig. lU^a)], such Fourier maps are rather 
featureless and only very low intensity is visible in the 
paths connecting nearest neighbor Nal-Na2 sites (note 
that as a result of the three inequivalent distances I, m h 
s the weak maxima are asymmetric). This suggests that 



the Na-ions are virtually static at these temperatures. 
At T ^ 200 K the inter-Na intensities slowly start to 
increase, in line with the anomaly of Nal <U22> mean- 
square displacement as well as previous NMR and /i+SR 
studies [lol, [HI- However, it is not until Ti is reached 
that a very strong increase in the residual scattering 
density is observed. As seen from Fig. lU^b), clear max- 
ima appear between the Nal-Na2 ions along the short 
5— paths, while along the long /—path the intensity re- 
mains weak. This is in excellent agreement with the 
picture presented above where Na-ion diffusion devel- 
ops through QID zigzag channels along the a-axis [see 
Fig. m^b)]. It should also be noted that there is a small 
intensity increase at (1/2, 0.85, 0.25) and equivalent po- 
sitions [labeled Na2' in Fig. lU^b)]. This Na-site is empty 
in the low and high temperature phases but becomes par- 
tially occupied in the orthorhombic phase due to the dis- 
placement of the Na2 position out of center of the triangle 
formed by the Nal sites. This displacement results in a 
small enlargement of the Na2' cavity, which explains the 
increase in the ab-plane area and unit cell volume for Ti 
> T > T2 mentioned in previous sections (see also [SM]). 
Rietveld refinements carried out including the Na2'-site 
indicate that it hosts about 2% of the Na+ ions between 
Ti and T2, which would correspond to the static part of 
the 'missing' Na-occupancy shown in Fig. [2fe-f). 

With further increasing temperature, another redistri- 
bution of the residual scattering intensity is observed at 
T2 [see Fig. lU^c)]. As the sample enters the hexagonal 
phase and all Nal-Na2 paths become equivalent (m), the 
previously closed channel along the b^/^— axis opens up 
and the Na-ion diffusion becomes isotropic within the 
ab-plane. Such QlD-to-2D evolution of the Na-diffusion 
mechanism is fully consistent with the above mentioned 
scenario suggested by the temperature dependence of the 
interatomic distances and anisotropic DW factors. Ergo, 
for this compound, structural transitions and Na dynam- 
ics are strongly interrelated. Although this behavior is 
found in other ionic conductors, the present case is par- 
ticular because Na diffusion not only starts along the 
shorter pathways: it also seems to trigger the structural 
transition which creates them, or viceversa. As demon- 
strated by our analysis, the activation of Na diffusion in 
a perfectly hexagonal lattice requires relatively high tem- 
peratures (T2 ^ 400 K). Hence, the occurrence of a very 
small structural distortion enables QID diffusion paths 
to open and Na-ion diffusion may set-in at much lower 
temperatures (Ti ^ 290 K). 

These findings are not only crucial for the understand- 
ing of the Nao.7Co02 compound, but also indicate that 
other materials from the same family may likely share a 
similar precursor to the 2D diffusive state. Consequently, 
it should be possible to change the onset temperature for 
diffusion by applying external perturbations {e.g. pres- 
sure) able to enhance the material's own intrinsic driv- 
ing force to create active diffusion channels. This offers 
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(a) T = 50 K « Ti (b) Ti < T = 350 K < T2 (c) T2 < T = 450 K 
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FIG. 4: (Color online) Fourier difference maps of the z = 0.25 Na planes at T = 50, 320 and 450 K showing the evolution of 
the residual scattering density in the paths connecting the Nal and Na2 sites (scale top/bottom: 0.25/0.0875 Fermi). For the 
three temperature regions there is a clear evolution from (a) Static Na ions, to (b) Quasi- ID ion diffusion along the ao-axis to 
finally (c) fully 2D ion diffusion. Thick solid white lines are overlayed together with the Na-ion positions in order to emphasize 
the diffusion paths. The shape of the Na ions has been adapted to the thermal ellipsoids previously shown in Fig. [3] and the 
different Nal-Na2 distances (I, m ,s) from Fig. [2{a) are also indicated. 



completely new and intriguing opportunities for an opti- 
mized tuning of functional materials in applications such 
a solid electrolytes or electrodes for novel rechargeable 
batteries. 

We finally address the impact of Na-ion diffusion on 
the C0O2 layers. As mentioned previously, no significant 
anomalies are observed in the Co-0 distances. However, 
the oxygen out-of-plane mean-square displacement 2x33 
undergoes a small but sharp decrease at Ti [see Fig.[3fe)]. 
This indicates that changes in Na landscape and the si- 
multaneous onset of QID Na diffusion indeed modify the 
width of the Co-0 distance distribution along the c-axis 
even if the mean value remains unchanged. Such behav- 
ior may arise from a modification of a static Na- vacancy 
order but more likely it is an effect related to the on- 
set of Na+ diffusion since the O atoms, closely bound to 
Co, do not follow easily the fast Na+ motion. Similar 
effects have been observed in magnetoresistive mangan- 
ites, where the crossover from adiabatic to non-adiabatic 
Mn-0 bond relaxation across the paramagnetic-insulator 
to ferromagnetic-metal transition also gives rise to a de- 
crease of the mean-square displacements [2o[. A dy- 
namic origin would also explain the dramatic changes 
of the ^^Co relaxation rate at Ti reported by NMR and 
/i+SR studies, whose unusual temperature dependence 
may simply reflect the state of the Na layers. 11, 2l| 

In summary, our data confirm the existence of struc- 
tural changes in Nao.7Co02 at Ti = 290 K and reveal 
another lattice anomaly at Ti = 400 K. A careful analy- 
sis of the crystal structure, the anisotropic Debye- Waller 
factors and the residual scattering density strongly sug- 
gest that both anomalies reflect a response of the lattice 
to the increasing Na mobility. The structural modifica- 
tions at Ti shorten 2/3 of the Nal-Na2 distances, creating 
QID pathways which allow for Na diffusion at temper- 
atures as low as 290 K. A full 2D diffusion is only ob- 
served above 400 K, coinciding with the stabilization of 



equal Na-Na distances in the hexagonal structure. The 
" melting" within the Na layers hence occurs in two steps 
and involves a crossover from QID to 2D Na diffusion. 
This clarifies two previously incompatible points of view, 
which associated the structural anomaly at Ti either to 
modifications of the static Na vacancy clustering or to 
the 2D-melting of the Na layers. It also shows that tai- 
lored ionic diffusion may be used to modify the crystal 
structure and vice versa. Such insight opens the door 
to the tuning of physical properties e.g. magnetism or 
electronic structure, and at the same time increases the 
possibilities for chemical design of novel materials with 
improved functional properties. 
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